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ABSTRACT 

Meteor showers on the Earth occur at well known times, and are associated with the decay 
of comets or other minor bodies whose orbital paths pass close to the Earth's trajectory. On the 
surface, determining the closest proximity of two orbital paths appears to be a computationally 
intensive procedure. This paper describes a simple geometric method for determining the prox- 
imity between the orbital paths of two bodies (i.e., a comet and a planet) in the solar system 
from the known ephemerides of the objects. The method is used to determine whether or not 
meteor showers on other planets in the solar system could be associated with any of 250 known 
comets. 

Subject headings: comets: general — meteors — solar system: general 



1. INTRODUCTION 

As they traverse their orbits about the Sun, 
comets slowly evaporate and fragment, leaving 
small bits of cometary debris along their orbital 
tracks. Some comet orbits intersect the Earth's 
path, and the planet sweeps up a portion of these 
particulates each year. Generally, these particles 
are drawn into the atmosphere, where they burn 
up at high altitudes, producing the yearly me- 
teor showers. A sample of the meteor showers 
expected on a regular basis for Earth-bound ob- 
servers is given in Table 1. A very detailed list 
of meteor streams encountered by the Earth has 
been composed based on ground-based observa- 
tions of amateur astronomers around the world 
(Jenniskens 1994). 

Given the large number of meteor showers seen 
on the Earth, it seems natural to ask about the 
possibility of meteor showers on other planets. It 
may be impractical for a sky-observer of the future 
to view meteor showers from some worlds: Mer- 
cury has no atmosphere, the clouds of Venus are 
so thick most meteors will likely burn up before 
a planetbound observer could see them, Jupiter 
has no solid surface to sit on while viewing the 
shower, and so forth. Never-the-less, predict- 



ing regular meteor showers on other worlds may 
be important for protecting explorers and space- 
craft from incoming particles, and could be use- 
ful for planning expeditions and experiments to 
collect cometary material. It has been suggested 
(Adolfsson et al. 1996) that future Mars landers 
may be able to detect meteors from the surface of 
Mars. 

A great deal of modern research has been 
devoted to analysis of the evolution of meteor 
streams in the solar system, particularly those 
that intersect the Earth's orbit (for example, de- 
tailed analyses of the evolution of the Quadran- 
tid stream can be found in (Hughes et al. 1979), 
(Williams et al. 1979) and (Murray et al. 1980); 
the Geminid stream is analyzed in (Fox et al. 1982) 
and (Jones 1985)). These analyses take into ac- 
count perturbations to the orbits of the parent 
bodies, as well as the subsequent evolution of the 
debris trail after the comet or minor body has 
continued on in its orbit. Over time, streams 
may wander into a planet's path causing new me- 
teor storms, or may wander out the planet's path 
quenching a shower which has been periodic for 
decades or centuries (e.g., (Murray et al. 1980) 
estimates that the Quadrantid shower will vanish 
by the year 2100). 
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To a first approximation, however, meteor 
showers will occur if the orbit of a planet and 
the orbit of a minor body intersect (or pass close 
to one another) . One way to determine if this oc- 
curs is to evolve the two orbits on a computer and 
watch for an intersection. A generalized method 
for finding the minimum separation between two 
orbits has been described by Murray, Hughes and 
Williams (1980), but the method reduces to a 
coupled set of equations for the orbital anomalies 
which requires numerical solution. Alternatively, 
the methods described in this paper approach the 
problem of determining the intersection of orbital 
planes in a completely analytical fashion, requir- 
ing only geometrical methods and matrix algebra. 

Section 2 describes the basic parameters and 
coordinate systems used to characterize orbits in 
this paper. Section 3 describes the rotations used 
to correctly orient two orbits with respect to each 
other, and applies the rotations to essential vec- 
tors needed for the analysis. Section 4 uses the 
rotated vectors to determine the intersection be- 
tween two orbital planes, and computes the dis- 
tance between the orbital paths when the planes 
intersect. Section 5 proposes a criteria for the ex- 
istence of a meteor shower based on the distance 
between the orbits at intersection. The "time" of 
showers meeting the criterion is determined. Sec- 
tion 6 applies the condition of Section 5 to 250 
known comets, summarizes the results, and dis- 
cusses the limitations of determining meteor show- 
ers using this method. 

Throughout this paper, SI (Systeme Interna- 
tionale) units are employed, except where the 
size of the units makes it convenient to work in 
standard units employed in astronomy (e.g., on 
large scales, astronomical units (AU) will be used, 
rather than meters). 

2. Describing Orbits 

As is well known, one of the great discoveries 
of Johannes Kepler was that the planets travel 
on elliptical paths, with the Sun at one focus of 
the ellipse (Kepler's First Law of Planetary Mo- 
tion, published in 1609). Since then, an enormous 
body of knowledge has been developed regarding 
the analysis of orbital motion (see, for example 
(Marion & Thornton 1988)), allowing the deter- 
mination of the position of virtually any object 



in the solar system at any moment in time. 

For the work presented here, a time dependent 
analysis of the orbital motion is not necessary 1 . 
The only information which is required is a knowl- 
edge of the trajectory of the orbit through space. 
The distance of the orbital path from the Sun may 
be written for elliptical orbits as 



where a is the semi-major axis of the orbit, e is the 
eccentricity, and 8 is the angle (the anomaly) be- 
tween the body and the axis defined by perihelion, 
as measured in the orbital plane. The perihelion 
distance for the object can be found from Eq. (1) 
by taking 9 = 0, yielding 

r p = o(l - e) . (2) 

The distance expressed in Eq. (1) describes the 
correct size and shape of an elliptical orbit for 
any object around the Sun, but more information 
is needed to correctly orient the orbit in three- 
dimensional space. This information is typically 
collected in the orbital ephemeris, which expresses 
orbital parameters with respect to the plane which 
is coincident with the orbital plane of the Earth 
(the ecliptic). This paper will use a (cartesian) 
reference coordinate system defined in the eclip- 
tic plane as shown in Figure 1. The +z axis is 
defined perpendicular to the ecliptic and in the 
right handed sense with respect to the Earth's or- 
bital motion (i.e., when viewed looking down the 
+z axis, the Earth's motion is counter-clockwise 
in the xy-plane) . The +x axis is defined along the 
direction of the Earth's perihelion. 

The orbital ephemeris of any body describes its 
orbit relative to the ecliptic plane, and locates the 
object along its orbital path as a function of time. 
For the problem of determining the possible inter- 
section of two orbital paths, only three elements 
of the full ephemeris for a body will be needed: 0, o 
(a modified longitude of the ascending node), t (in- 
clination), and lo (argument of perihelion). Each 
of these parameters is described below, and shown 
in Figure 2. 

1 We are interested in knowing the proximity of two orbits 
when they cross. An interesting (but ultimately more dif- 
ficult) question to address is whether two bodies might ac- 
tually collide because their orbits intersect. 
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The longitude of the ascending node, O, is 
the angle in the ecliptic plane between the vernal 
equinox (the first point of Aries) and the point at 
which the orbit crosses the ecliptic towards the +z 
direction ("northward" across the ecliptic). The 
parameter, Q a , used in this paper, is an offset 
longitude measured from the perihelion of Earth, 
rather than the first point of Aries (see Figure 3). 
The inclination, i, is the angle between the nor- 
mal vector of the orbit and the normal vector of 
the ecliptic. Lastly, the argument of perihelion, u>, 
is the angle between the position of the body as 
it crosses the ascending node and the position at 
perihelion, as measured in the orbital plane of the 
body. 

In addition to these three angles, it will be use- 
ful to define two vectors for each orbit of interest: 
h, the unit normal vector to the plane of the or- 
bit, and f p , the vector pointing to perihelion in 
the plane of the orbit. 

3. Rotations for Orbital Orientation 

In order to correctly orient an orbit with respect 
to the ecliptic, assume (initially) that the orbit of 
interest is in the plane of the ecliptic, with the 
perihelion of the orbit aligned along the +x axis 
(i.e., the orbit is co-aligned with the Earth's or- 
bit). A series of three rotations, based on the an- 
gles {u>, 6, fl} from the orbital ephemeris will pro- 
duce the correct orientation. The first rotation will 
set the value of the ascending node with respect 
to perihelion, the second rotation will set the incli- 
nation to the ecliptic, and the third rotation will 
move the ascending node to the correct location 
in the ecliptic plane. 

A useful method for describing rotations is in 
terms of matrices. While it is possible to construct 
a rotation matrix for rotations about a general 
axis, it is more convenient to conduct rotations 
about the coordinate axes shown in Figure 1. The 
matrices describing rotations about the x-,y-, and 
z-axes will be denoted M x (<j>), M y (£_), and M z (ip), 
respectively. 

To demonstrate the rotations needed to orient 
the orbit, consider a general vector, A, which is 
rigidly attached to the orbital plane, maintaining 
its orientation as the plane is rotated. 

The first rotation locates the ascending node 
with respect to perihelion; the rotation depends on 



the value of the argument of perihelion, u>. This is 
done by rotating around the z-axis by tjj = u). In 
terms of rotating a general vector A, this can be 
written 

A 1 =M z {uj)A. (3) 

When this operation is applied to the orbit, the 
ascending node will be located on the +x axis. 

The orbit is inclined around an axis which 
passes through the ascending node and through 
the Sun (at one focus of the orbit). Since the first 
rotation placed the ascending node on the +x axis, 
and the Sun lies at the origin of coordinates, a ro- 
tation around the x-axis by the inclination angle, 
<fi = l, will correctly incline the orbit. In terms of 
the vector A\ (resulting from Eq. (3)), this yields 



A 2 = M x (l)A 1 



(4) 



Before the final rotation, it will be convenient 
to offset the longitude of the ascending node such 
that it is measured from the perihelion of the 
Earth, rather than the vernal equinox (this makes 
the x-axis the origin for measuring the longitude 
of the ascending node). The angle between the 
vernal equinox and perihelion of Earth is simply 
the argument of perihelion for Earth, w®, giving 
(see Figure 3) 



f2 = 2-7T — u>® + Q, . 



(5) 



After the second rotation, the ascending node 
is still located on the +x axis. Rotation about 
the z-axis by the offset longitude, ip = £1 , will 
rotate the longitude of the ascending node to its 
correct location in the ecliptic plane. In terms of 
the vector A 2 (resulting from Eq. (4)), this yields 



A 3 = M Z (Q )A 2 



(6) 



The vector A3 (which is rigidly attached to the or- 
bit) is correctly oriented with respect to the eclip- 
tic. 

The two vectors which will be of use later are 
the unit normal vector to the orbit, n, and the 
perihelion vector, r p . When the orbital plane is 
co-aligned with the Earth's (before any rotations 
have been performed) , these vectors have the form 



'v 





(7) 
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The rotation operations described by Eqs. (3), (4), 
and (6) must be applied to these vectors so they 
correctly describe the orbit with respect to the 
ecliptic. Conducting the rotation procedure yields 



n = 



sin i sin Q Q 
- sin i cos fl 

COS L 



(8) 



and 



cos ui cos il a — sin ui cos i sin £1 
cos u> sin o + sin u cos t cos fl a 
sin w sin/, 



(9) 



4. Intersection of Orbits 



The procedure described in Section 3 will cor- 
rectly orient any orbit with respect to the ecliptic. 
One could take any planet's cphemeris (e.g., from 
the ephemerides given in Table 2) and construct 
the normal vector n and perihelion vector f p in ac- 
cordance with Eqs. (8) and (9) 2 . Similar vectors 
could be generated for cometary ephemerides. 

The real question of interest is not how the or- 
bital planes of planets and comets are related to 
the ecliptic, but rather how they are oriented with 
respect to each other, and in particular where they 
intersect. The line defining the intersection of the 
orbital planes can be used to determine whether 
or not the orbits actually intersect. 

Hereafter, assume that vectors related to a 
comet's orbit will bear the subscript 'c' and vectors 
related to a planet's orbit will bear the subscript 
'+'. Further, suppose the components of the nor- 
mal vector for a comet's orbit are n c = (a,b,c), 
and the components of the normal vector of a 
planet's orbit are h + = (e,f,g). Both orbital 
planes automatically share one point in common: 
the origin, which lies at the focus of each orbital 
ellipse. Given this point and the two vectors n c 
and n + , the equations describing the two orbital 
planes are 



CometPlane 
PlanctPlane 



ax + by + cz = 
ex + fy + gz = 



(10) 



2 To ease the notation, we will drop the primed notation for 
rotated vectors from here on. It will be understood that the 
normal vectors and perihelion vectors have been correctly 
oriented with respect to the ecliptic. 



The intersection of the two planes is a line 
which is the common solution of the two expres- 
sions in Eq. (10). Using determinants, the com- 
mon solution to these equations is found to be 



(11) 



X 




-y 




z 


b c 




a c 




a b 


f g 




e g 




e f 



where k is an arbitrary constant. The solutions 
{x,y,z} of Eq. (11) will be points along the line 
of intersection. It is useful to use these values to 
define a new vector, A, called the 'node vector.' It 
points along the line of nodes (the intersection of 
the two planes), and has components 



bg-cf 

ce — ag 
af — be 



(12) 



To determine if the orbital paths intersect, one 
must know the radii of the orbits along the line of 
nodes. An orbital radius may be determined from 
Eq. (1) if the value of the anomaly, 9, is known. In 
terms of two orbits inclined with respect to each 
other, the angles of interest will be the angle be- 
tween the perihelion vector for each orbit, f p , and 
the node vector, A. For each orbit, the angle is 
defined in terms of the dot product of the two 
vectors, yielding 



cos 6* 



f P ■ A 



K\ ■ A 



(13) 



The orbits have two opportunities to intersect: at 
the ascending node, and at the descending node. 
Eq. (13) gives the angle at a single node. To obtain 
the value of the anomaly at the other node, dot the 
perihelion vector, f p , into the negative of the node 
vector, —A. 

Once the anomaly is known, the distance be- 
tween the orbital paths when the planes intersect 
is simply 

A = |r+-r c |, (14) 

where r + and r c are computed using Eq. (1) with 
the anomaly defined by Eq. (13) and the appro- 
priate orbital parameters derived from tabulated 
ephemerides. 
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5. Is there a meteor shower? 

The occurrence of a meteor shower associated 
with a particular comet will depend on the value 
of the separation between the orbital paths, A. A 
variety of proximity criteria could be developed, 
depending on specific considerations one would 
like to make regarding the likclyhood that parti- 
cles from a given stream might reach a planetary 
body (such as the proximity of the mean orbits 
of a comet and planet, or the spread of cometary 
particles around the orbit of the parent comet). 
The basic criteria may be expressed such that 



A < f 



(15) 



where r is the minimum separation one considers 
likely for a shower to occur. 

In this paper, the criteria for an orbit inter- 
section causing a meteor shower will be scaled to 
a region where the gravitational potential of the 
planet dominates over the gravitational potential 
of the Sun. The criteria will be 



A < kRi , 



(16) 



where Ri is the "Roche-lobe radius" (defined as 
the radius of a sphere which has the same volume 
as the planet's Roche lobe) and n is an arbitrary 
adjustable scaling factor. The Roche-lobe radius 
can be approximated by 



R, -0.52 -a 



m+ 



n 0.44 



M Q + m + 



(17) 



where m + and Mq are the mass of the planet 
and the Sun, and a is the semi-major axis of the 
planet's orbit (Iben & Tutukov 1984). 

Once an intersection (in the sense of Eq. (16)) 
has been found, one would like to know when the 
associated meteor shower might occur, such that 
it is possible to mount an observational effort to 
detect the shower. A good marker for a plane- 
tary encounter with the meteor stream is the value 
of the planetary anomaly along the line of nodes, 
given by Eq. (13). 

Given a particular value of the anomaly, 9, it 
is possible to write a closed- form expression t(9) 
for the time at which the planet will arrive at that 
point in its orbit. The areal velocity may be writ- 
ten 

dA nab 

~dt = —' (18) 



where a and b are the semi-major and -minor axes 
respectively, and r is the period of the orbit. Sep- 
arating Eq. (18) and writing the area element as 
dA = (l/2)r 2 d9 allows one to write an integral for 
the time t as 



t{6) 



(—)- 

KirabJ 2 



r 2 d9' . 



(19) 



Using the shape equation (Eq. (1)) to write r = 
r{9), and using b = o(l — e 2 ) 1 / 2 to express the 
semi-minor axis, Eq. (19) can be integrated to give 



t{9) = 



2tt 



2 tan 




ey/1 — e 2 sin0 



1 + e cos 9 
(20) 

which is the time it takes the planet to pass from 
perihelion (9 = 0) to the anomaly value 9. 

6. Results & Discussion 

As an example of the methods presented 
here, a search for comet-planet orbital intersec- 
tions in the solar system was carried out us- 
ing the planetary ephemerides shown in Table 2 
(Standish et al. 1992), and the comet ephemerides 
provided in the Jet Propulsion Laboratory's 
DASTCOM (Database of ASTeroids and COMets) 
(JPL DASTCOM 1999). The DASTCOM is a col- 
lection of orbital parameters and physical charac- 
teristics for the numbered asteroids, unnumbered 
asteroids, and periodic comets, used for analyses of 
solar system dynamics. The DASTCOM includes 
all known comets with periods less than 200 years 
(periodic comets), and long period comets which 
have made a return after 1995 (there are cur- 
rently 73 long period comets in the DASTCOM, 
and 208 comets with periods less than 200 years) . 
The orbital elements given in the DASTCOM are 
osculating elements, computed from observations 
during the comet's most recent apparition. 

The results of the search for comet-planet or- 
bital intersections are listed in Table 3, which spec- 
ified an encounter distance of 



A < 5Ri . 



(21) 



In all, 128 possible showers were detected: 3 at 
Earth, 1 at Mars, 106 at Jupiter 3 , 17 at Saturn, 



3 In fact, the results of Table 3 show that Jupiter's orbit 
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and 1 at Uranus. If one reduces the encounter 
distance to A < lRi, only 32 possible showers 
are detected (shown at the top of Table 3): 1 at 
Earth, 28 at Jupiter, 2 at Saturn, and 1 at Uranus. 
If one allows the encounter distance to expand to 
A < 10i?;, 188 possible showers are detected (data 
not shown in Table 3): 4 at Earth, 5 at Mars, 148 
at Jupiter, 24 at Saturn, 6 at Uranus, and 1 at 
Neptune. 

Comets with orbital scales smaller than the so- 
lar system ('short period comets') have evolved 
largely under the influence of perturbations due 
to Jupiter (the mass of Jupiter is greater than 
the mass of the other planets combined), giving a 
large population of comets which cross Jupiter's 
orbit. The search for the origin of these "Jo- 
vian family comets" has been a matter of much 
numerical simulation and debate (see, for exam- 
ple, (Quinn et al. 1990)). The disproportionately 
large number of showers detected for Jupiter can 
be attributed to this feature of the comet popula- 
tion. 

A good check of the procedure described in 
this paper is to consider the predicted showers at 
Earth. In particular, the method outlined in this 
work predicts two meteor streams which can be 
identified with known showers. The first is the 
stream from Comet Tempel-Tuttle, occurring at 
t <~ 318 d. This stream can be identified with the 
Leonid meteor shower (known to be a stream from 
Tempel-Tuttle), which occurs in mid-November 
each year. The second is a stream from Comet 
Swift-Tuttle, occurring at t <~ 221 d. This stream 
can be identified with the Perseid meteor shower 
(known to be a stream from Swift-Tuttle), which 
occurs in mid- August each year. 

The possible showers computed here have all as- 
sumed that the orbits of the comets are static and 
do not precess. Further, it is assumed that the me- 
teor streams remain attached to those static orbits 
without wandering under the influence of gravita- 
tional perturbations in the solar system. In addi- 
tion, the influences of 'local' bodies around each 
planet (e.g., Earth's moon, or the Galilean satel- 
lites around Jupiter) have been ignored. Never- 

intcrsccts the path of comet P/Spahr (1998 U4) twice: one 
intersection at a separation of A ~ 1.6-R;, and a second 
intersection (at the other node) with a separation of A ~ 
3.3R ; . 



thc-less, the method provides a useful way for de- 
termining the possibility that a given planet will 
encounter a meteor stream from minor bodies in 
the solar system. 

I would like to thank M. B. Larson for com- 
ments and suggestions, and E. M. Standish who 
provided helpful discussions regarding planetary 
ephemerides. D. K. Ycomans and P. Chodas 
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database. I would also like to acknowledge the 
hospitality of the Solar System Dynamics Group 
at the Jet Propulsion Laboratory during the time 
this work was completed. This work was sup- 
ported in part by NASA Cooperative Agreement 
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Fig. 1. — The reference coordinate system in the 
ecliptic plane. The z-axis is defined in the right- 
handed sense with respect to the Earth's motion, 
and the x-axis points towards the perihelion of the 
Earth. 



Fig. 2. — The three essential angles for correctly 
orienting orbits in three dimensional space are (a) 
o , the (modified) longitude of the ascending node; 
(b) l, the inclination; and (c) ui, the argument of 
perihelion. 



Fig. 3. — The modified longitude of the ascending 
node, o , defined in terms of the Earth's argument 
of perihelion, , and the conventional value of f2 
for the orbit, v indicates the first point of Aries. 



Fig. 4. — The intersection of two orbits, showing 
the essential quantities for defining the occurrence 
of a meteor shower: the separation of the orbits at 
crossing, A and the node vector, A, which defines 
the intersection of the two orbits and is used to 
determine the time of the meteor shower. 
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Table 1 

Some yearly meteor showers seen from Earth. 



Shower Name 


Date 


Quadrantids 


early January 


Lyrids 


mid April 


T] Aquarids 


early May 


5 Aquarids 


late July 


Perseids 


mid August 


Orionids 


mid October 


Leonids 


mid November 


Geminids 


mid December 



Table 2 

The mean ephemerides for the planets of the solar system (epoch J2000). 



Planet 


a 


e 


i 


n 






(AU) 




(°) 


(°) 


(°) 


Mercury 


0.38709893 


0.20563069 


7.00487 


48.33167 


29.12478 


Venus 


0.72333199 


0.00677323 


3.39471 


76.68069 


54.85229 


Earth 


1.00000011 


0.01671022 


0.00005 


-11.26064 


114.20783 


Mars 


1.52366231 


0.09341233 


1.85061 


49.57854 


286.4623 


Jupiter 


5.20336301 


0.04839266 


1.3053 


100.55615 


-85.8023 


Saturn 


9.53707032 


0.0541506 


2.48446 


113.71504 


-21.2831 


Uranus 


19.19126393 


0.04716771 


0.76986 


74.22988 


96.73436 


Neptune 


30.06896348 


0.00858587 


1.76917 


131.72169 


-86.75034 


Pluto 


39.348168677 


0.24880766 


17.14175 


110.30347 


113.76329 
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Table 3 

The results of a meteor shower search using the comets in the JPL DASTCOM database. 



Planet 


Comet 


(A/R ; ) a 


8, degrees 


t{6),d 


Earth 


109P/Swift-Tuttle 


0.5043 


216.49717 


220.81696 


Jupiter 


P/LONEOS-Tucker (1998 QP54) 


0.02759 


143.72224 


1688.24555 


Jupiter 


117P/Helin- Roman- Alu 1 


0.05761 


54.79245 


605.7507 


Jupiter 


43P /Wolf-Harrington 


0.10689 


241.59541 


2965.94337 


Jupiter 


P/Hergenrother (1998 W2) 


0.13008 


158.6528 


1883.37507 


Jupiter 


C/Hale-Bopp (1995 Ol) 


0.15907 


267.71691 


3287.20798 


Jupiter 


78P/Gehrels 2 


0.2102 


206.19238 


2510.80676 


Jupiter 


75P/Kohoutek 


0.21718 


256.88227 


3155.62992 


Jupiter 


P/Spahr (1998 Wl) 


0.24054 


267.30471 


3282.24601 


Jupiter 


124P/Mrkos 


0.32574 


344.80757 


4164.72657 


Jupiter 


14P/Wolf 


0.32874 


191.77241 


2321.02269 


Jupiter 


53P/Van Biesbroeck 


0.37063 


143.86166 


1690.05544 


Jupiter 


59P/Kcarns-Kwee 


0.39771 


294.48157 


3602.23635 


Jupiter 


91P/Russell 3 


0.41994 


56.45964 


624.68681 


Jupiter 


76P / Wcst-Kohoutek-Ikemura 


0.44039 


248.71583 


3054.87529 


Jupiter 


26P/Grigg-Skjellerup 


0.45251 


21.59907 


236.07478 


Jupiter 


132P/Helin- Roman- Alu 2 


0.50932 


176.74286 


2122.18424 


Jupiter 


P/Kushida (1994 Al) 


0.54051 


239.05142 


2933.93263 


Jupiter 


16P/Brooks 2 


0.57961 


175.77613 


2109.39006 


Jupiter 


83P/Russell 1 


0.59462 


34.84717 


382.18644 


Jupiter 


D/Kowal-Mrkos (1984 HI) 


0.61747 


57.35711 


634.89823 


Jupiter 


135P/Shoemaker-Levy 8 


0.68357 


28.92759 


316.71554 


Jupiter 


139P /Vaisala-Otcrma 


0.72789 


241.18336 


2960.76693 


Jupiter 


86P/Wild 3 


0.73742 


55.46464 


613.3804 


Jupiter 


104P/Kowal 2 


0.81618 


233.8782 


2868.46755 


Jupiter 


P/LINEAR-Mueller (1998 SI) 


0.87419 


157.76011 


1871.64355 


Jupiter 


P/Shoemaker-Levy 6 (1991 VI) 


0.87826 


199.12624 


2418.00496 


Jupiter 


18P/Perrine-Mrkos 


0.94922 


228.33725 


2797.82077 


Jupiter 


85P/Boethin 


0.97138 


164.48654 


1960.19405 


Saturn 


P/Jager (1998 U3) 


0.06673 


210.03052 


6366.0197 


Saturn 


126P/IRAS 


0.53161 


83.17079 


2299.91062 


Uranus 


C/Li (1999 El) 


0.84435 


137.45475 


11360.58497 


Earth 


55P/Tcmpcl-Tuttlc 


4.15373 


312.31161 


318.28395 


Earth 


26P/Grigg-Skjellerup 


4.518 


110.36191 


110.13972 


Mars 


C/LINEAR (1998 U5) 


1.98779 


90.14579 


151.61158 


Jupiter 


47P / Ashbrook- Jackson 


1.00145 


162.12665 


1929.08903 


Jupiter 


15P/Finlay 


1.00753 


186.90971 


2256.74179 


Jupiter 


97P/Mctcalf-Brcwmgton 


1.09545 


175.71325 


2108.55792 


Jupiter 


81P/Wild 2 


1.10168 


320.59158 


3897.71867 


Jupiter 


121P/Shoemaker-Holt 2 


1.10447 


264.90153 


3253.24816 


Jupiter 


54P/dc Vico-Swift 


1.1133 


152.93248 


1808.3284 


Jupiter 


56P/ Slaughtcr-Burnham 


1.1211 


143.93726 


1691.03684 


Jupiter 


P/Korlevic-Juric (1999 DN3) 


1.15484 


347.38897 


4192.9561 


Jupiter 


P/Mucllcr 4 (1992 G3) 


1.16212 


312.27537 


3804.66421 


Jupiter 


52P/Harrington-Abell 


1.20057 


316.89669 


3856.48162 


Jupiter 


69P/Taylor 


1.2392 


274.64237 


3370.06049 


Jupiter 


P/Larsen (1997 VI) 


1.269 


224.13556 


2743.90598 


Jupiter 


46P/Wirtanen 


1.31801 


245.23688 


3011.54806 


Jupiter 


lOOP/Hartley 1 


1.32008 


21.70502 


237.23779 


Jupiter 


87P/Bus 


1.34382 


15.87435 


173.3314 


Jupiter 


C/Fcrris (1999 K2) 


1.3836 


285.50894 


3498.14812 


Jupiter 


77P/Longmore 


1.43731 


358.02608 


4309.06758 


Jupiter 


C/Mueller (1997 Jl) 


1.44651 


262.26797 


3221.3349 


Jupiter 


P/Hartley-IRAS (1983 VI) 


1.51103 


166.77415 


1990.37975 
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Table 3 — Continued 



Planet 


Comet 


( A /Ri) 


0, degrees 


t{0), a 


Jupiter 


119P /Parker-Hartley 


1.53129 


236.51823 


2901.93741 


Jupiter 


114P /Wiseman-Skiff 


1.54424 


256.92325 


3156.13216 


Jupiter 


102P/Shoemaker 1 


1.58558 


143.00259 


1678.90784 


Jupiter 


P/Spahr (1998 U4) 


1.61195 


169.07697 


2020.79437 


Jupiter 


33P/Danicl 


1.64668 


229.96736 


2818.65939 


Jupiter 


62P/Tsuchinshan 1 


1.6969 


261.53064 


3212.37462 


Jupiter 


70P/Kojima 


1.78536 


288.99826 


3538.80021 


Jupiter 


60P/Tsuchinshan 2 


1.84257 


272.21092 


3341.08206 


Jupiter 


4P /Faye 


1.91557 


192.45922 


2330.09358 


Jupiter 


67P /Churyumov-Gerasimcnko 


1.95263 


207.28864 


2525.16175 


Jupiter 


6P /d'Arrest 


2.01111 


126.65145 


1468.85806 


Jupiter 


36P /Whipple 


2.02544 


175.2223 


2102.06112 


Jupiter 


C/LINEAR (1998 Ul) 


2.08822 


12.83484 


140.08603 


Jupiter 


C/Spacewatch (1997 BA6) 


2.15772 


302.67334 


3696.04414 


Jupiter 


P/Levy (1991 L3) 


2.21811 


131.97724 


1536.80051 


Jupiter 


116P /Wild 4 


2 27558 


347 04647 


4189.21226 


Jupiter 


P/Shoemaker-Levy 1 (1990 VI) 


2.39831 


215.04783 


2626.37534 


Jupiter 


9P /Tempcl 1 


2.41105 


50.10947 


552.78419 


Jupiter 


C/LINEAR (1999 H3) 


2.59549 


138.47987 


1620.39097 


Jupiter 


31P/Schwassmann-Wachmann 2 


2.63332 


303.86438 


3709.59249 


Jupiter 


P/LINEAR (1999 J5) 


2.69947 


98.43266 


1117.77483 


Jupiter 


C/LINEAR (1998 W3) 


2.71291 


108.75604 


1244.38557 


Jupiter 


21P/ Giacobini-Zinncr 


2.73222 


182.72852 


2201.41502 


Jupiter 


40P/Vaisala 1 


2.84407 


304.29294 


3714.46212 


Jupiter 


108P /Ciffreo 


2.88171 


214.58058 


2620.30119 


Jupiter 


136P/Mueller 3 


2.92721 


128.57415 


1493.32993 


Jupiter 


42P/Neujmin 3 


2.93144 


153.24921 


1812.47528 


Jupiter 


7P /Pons-Winnecke 


3.13747 


78.25359 


876.539 


Jupiter 


103P /Hartley 2 


3.19169 


209.67222 


2556.32919 


Jupiter 


128P /Shoemaker-Holt 1-B 


3.22265 


213.44774 


2605.56299 


Jupiter 


D/van Houten (1960 SI) 


3.25993 


336.58556 


4074.58242 


Jupiter 


P/Kushida-Muramatsu (1993 XI) 


3.26773 


250.65642 


3078.93931 


Jupiter 


C/Zhu-Balam (1997 LI) 


3.26781 


218.82993 


2675.43314 


Jupiter 


P/LINEAR (1998 VS24) 


3.26871 


158.78542 


1885.11843 


Jupiter 


P/Spahr (1998 U4) 


3.29681 


349.07697 


4211.40146 


Jupiter 


61P /Sliajn-Schaldach 


3.3833 


164.21702 


1956.63978 


Jupiter 


65P /Gunn 


3.43747 


49.5381 


546.34347 


Jupiter 


120P/Mueller 1 


3.59615 


161.60426 


1922.20873 


Jupiter 


129P/Shoemaker-Levy 3 


3.6607 


284.2452 


3483.36888 


Jupiter 


30P/Reinmuth 1 


3.66617 


288.52425 


3533.29093 


Jupiter 


22P/Kopff 


3.67008 


113.53434 


1303.71706 


Jupiter 


112P/Urata-Niijima 


3.69814 


194.43253 


2356.14237 


Jupiter 


HOP/Hartley 3 


3.70194 


272.28798 


3342.00224 


Jupiter 


49P / Arend-Rigaux 


3.78256 


288.50492 


3533.06618 


Jupiter 


P/Lagerkvist (1996 R2) 


3.9786 


175.43706 


2104.90302 


Jupiter 


19P/Borrelly 


4.00335 


239.67933 


2941.845 


Jupiter 


17P/Holmcs 


4.00357 


130.61884 


1519.42495 


Jupiter 


137P/Shoemaker-Levy 2 


4.01952 


229.51888 


2812.93064 


Jupiter 


131P/Mueller 2 


4.05027 


208.10376 


2535.82718 


Jupiter 


P/Jedicke (1995 Al) 


4.07511 


102.29178 


1164.85082 


Jupiter 


D/Tritton (1978 C2) 


4.234 


282.865 


3467.19321 


Jupiter 


48P/ Johnson 


4.31823 


104.30891 


1189.57815 


Jupiter 


106P/Schuster 


4.41102 


212.99144 


2599.62197 


Jupiter 


C/LINEAR (1998 M5) 


4.67076 


318.4884 


3874.26616 
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Table 3 — Continued 



Planet 


Comet 


(A/R t r 


8, degrees 


t(9),d 


Jupiter 


P/LONEOS (1999 R028) 


4.84285 


141.12541 


1654.58458 


Jupiter 


98P /Takamizawa 


4.88896 


113.76163 


1306.55052 


Jupiter 


C/Spaccwatch (1997 P2) 


4.89278 


286.39551 


3508.49848 


Jupiter 


73P/Schwassmann-Wachmann 3 


4.9158 


51.54895 


569.03117 


Jupiter 


P/Helin-Lawrencc (1993 K2) 


4.99431 


76.00361 


850.14018 


Saturn 


C/Catalina (1999 Fl) 


1.01329 


107.67064 


3035.62842 


Saturn 


P/Gehrels (1997 CI) 


1.27361 


164.60138 


4862.59895 


Saturn 


P/Hermann (1999 Dl) 


1.35745 


252.08077 


7703.6729 


Saturn 


C/LINEAR (1998 Ql) 


1.36923 


70.34971 


1928.09742 


Saturn 


P/Shoemaker 4 (1994 J3) 


2.45512 


358.40089 


10704.16966 


Saturn 


P/Montani (1997 Gl) 


2.89611 


185.28975 


5549.17747 


Saturn 


P/Helin (1987 Q3) 


3.16833 


76.81205 


2114.42428 


Saturn 


63P/Wild 1 


3.2706 


260.23561 


7952.44828 


Saturn 


140P/Bowell-Skiff 


3.30576 


231.84877 


7070.65503 


Saturn 


D/Bradfield 1 (1984 Al) 


3.57748 


262.73062 


8027.80737 


Saturn 


134P/Kowal- Vavrova 


3.64885 


319.9063 


9665.71078 


Saturn 


C/Spacewatch (1997 BA6) 


3.95068 


44.94144 


1214.46702 


Saturn 


C/LINEAR (1999 N4) 


3.95944 


78.38428 


2160.08505 


Saturn 


C/LINEAR (1999 H3) 


4.19679 


241.0887 


7362.46033 


Saturn 


P/Lagcrkvist-Carscnty (1997 T3) 


4.90066 


120.32592 


3428.86878 



a The first 32 entries (above the line) are for intersections having A < Ri. All other 
encounters are for A < 5Ri. 



11 



" +z 



ascending 
node 



